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ABSTRACT 


An innovative thin film deposition technique : low 
energy ion beam deposition , is further studied. The 
theoretical deposition processes of this method are 
discussed. Double layer thin films of Mg and Pb on carbon 
substrates are produced by this method and compared with 
Similar films grown by thermal evaporation. Their 
characteristics are investigated by Rutherford 
backscattering analysis and Scanning Electron Microscopy. 
The RBS technique is expanded to derive a profile of 
composition ratio versus depth within the deposits and 
applied to a depth analysis of the double layer films. The 
present study of film characteristics shows that a high 
level of mass selectivity is observed from the low energy 
ion beam deposited Mg films. The oxidation of Mg films is 
not limited to the film surface : oxygen is seen to be 
distributed into the Mg region rather than staying on the 
surface. The preferential migration of Pb into Mg is 
observed and the diffusion coefficient of Pb into Mg at 300 
C is calculated as 1.236 x MO ccuyseee The interdiffusion 
region between the Mg and Pb film layer ranges from 1000 to 
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3000 A at room temperature, 
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CHARTERS 4 


INTRODUCTION 


The technology and understanding of thin films has 
advanced tremendously in recent years because of an 
industrial demand for microelectronic devices and optical 
coatings. In 1852 Bunsen and Grove (*1) obtained the first 
thin films by chemical reaction and glow-discharge 
sputtering techniques. Faraday (*1) also produced metal 
films by thermal evaporation in 1857. Since these initial 
efforts the techniques of thermal evaporation and glow 
discharge sputtering have been most frequently used for thin 
film production. With the advent of high power laser 
technology the need for high damage threshold material is 
apparent. 

The damage to the optical interface is considered the limit 
of laser performance and the properties of the optical 
interface is limited by the film inhomogeneities, non linear 
susceptibility or weak absorption (*2). Various kinds of 
films have been made for laser optics. Most of these films 
are made by deposition of film material by evaporation or by 
sputtering onto a suitable substrate material. In both 
methods the substrate is within a high pressure environment 


as is the source of the coating material, a procedure which 
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increases the risk of contamination . The production of thin 
films by using a low energy, mass analyzed metallic ion bean 
under ultra high vacuum conditions was successfully 
demonstrated by Amano et. el. (*3). The characteristics of 
Single layér films of lead and magnesium deposited on the 
Carbon substrate by a low energy ion beam technique was 
investigated by Amano et. al. (*4). 

In the low energy ion beam deposition method, the coating 
material itself is converted into an ion beam which is mass 
and energy analyzed and then transported to a differentially 
pumped chamber in which the substrate resides under ultra 
high vacuum conditions. Film growth proceeds when the ion 
beam is retarded to an energy low enough to minimize 
substrate or film damage but high enough to produce a good 
bond to the substrate. 

The controllability of the ionized beam contributes 
favorably to the deposited film characteristics. A detailed 
comparison is outlined in Table 1-1, 

The purpose of this thesis is to investigate 
characteristics of low energy ion beam deposited double 
layer films of lead and magnesium and to compare them with 
those made by a thermal evaporation method. 

A Rutherford backscattering analysis method is chiefly 
employed to analyze the films. Film surface structure is 
determined by a scanning electron microscope and the 
chemical composition of magnesium films is confirmed by 


ESCA (Electron Spectroscopy for Chemical Analysis). 
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A computer program is written to translate the 
RBS(Rutherford BackScattering) spectra into physical depth 
and atomic concentration profiles. 

The following example explains the convention used in 
this thesis to identify the film material in order of 
deposition and the method of deposition. 

A number in parentheses following a letter identifies that 
material having been deposited by the IBD technique at the 
aera (in eV) denoted by the number. The letter T in 
parentheses following a letter identifies a thermally 
evaporated film. The order of the letters after C (the 


carbon substrate) identifies the order of deposition. 


ex) CM (48)P(T) 


SS — ee ee oe ee 


carbon Mg Pb 


substrate 


CM (48)P(T) identifies a carbon substrate with a magnesium 
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film ion beam deposited at 48 eV energy, follwed by a 


thermally evaporated lead layer. 
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CHAPTER 2 


THE LOW ENERGY ION BEAM DEPOSITION PROCESS 


Ze 1 Introduction 


Starting from the energetic ions moving towards a 


target until the film forms on the conductive target, the 


following processes are considered to be happening. The 


following is just a summary of the processes and they will 


be discussed in more detail in later sections. 


a 


The ions are neutralized near the conductive 
target surface before the kinetic collison process 
can occur. 

As the neutralized atom moves closer to the target 
surface, the kinetic energy of the neutralized 
atom is transferred to the target by elastic 
collisions because the initial energy of the 
incident ion is low (*1). 

Losing kinetic energy, the adatom will be 
accommodated to the surface causing nuclei to 
form; the nuclei grow in three dimensions forming 
islands and the islands grow until they touch. The 
neighboring islands coalesce, forming larger 


islands. As the islands continue to grow, the 
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tendency of maintaining roundness after 
coalescence decreases. The islands become 
elongated and join to form continuous networks 
leaving long and narrow channels between them. As 
the deposition proceeds, the channels are filled 
and a continuous film is formed. 

d ; During and after the deposition process, diffusion 
of the film material occurs in three dimensions. 

e ; When the kinetic energy of the ion is higher than 
the threshold energy for displacement of the film 
material, self-sputtering adversely affects the 


film growth. 


2.2 Neutralization of ions on the metal surface and emission 


of secondary electrons. 


In considering the interaction between an ion anda 
metal atom when an ion is approaching the metal surface, it 
is convenient to distinguish between the kinetic energy and 
potential energy of the ion. The kinetic energy is provided 
by the acceleration potential, usually applied between the 
ion source and target, while the potential energy is given 
by the ionization process. 

For low energy incoming ions interacting with target 
material, the transfer of kinetic energy from the former to 


atoms of the latter can be adequately described using hard 


ebnnte: 


‘ayona 2 tabs 


ot 3 i Jae 23 Be sweeie 


a23% 


mes Aner 
mld “* ‘foe 8 
a ak sods seen weep 


i 
boll se 976 aflegasdo | sit n2bGs 9834 


Pay 


428 5634 foitr 


«bewxed 2 a@it> 


_ 


7a 


§ >» 
sopeveety’ 


. 


2094 y Oo 7 


-erotanesid sea¢s¢ ai exzvg00 1hra 
7 i = 7 
tedd tetiohd sk pot ed?) to yotend ae A oft A948 * 
om 
ali €A7% +o tazasoshqeeh 1 eit 6 foie sida) siz 
"a | —— Per , 
+ e2997Ts yiectov is a rtesroqge-2tes  ~lbi3Siam sl | 
ae ~ 
: iW Oy R afst 7 — 
a 
7 = . 

_ 7 - ~s 
soiesazige ba si1bpe fetem add go enot Io aobstetiarznet ¢.f 
Diet 2 fed ait 3 rst 

| ; ’ r = > ~—— > 
= eaallQiI7 39812 ¥° 5586 2 TO 
i 7 _ 
i y Teor wu: y ; 7 
; ny ‘ won 4 ; 7 _ on 
7 5 ie, Af. 
has mob a5 teeussd. aoktsaaes, palzebée noo at 
>) ; i : 7 : 
: ,sbetuwe Letom si} pirntdssotag i ts cade cots Leaves 
°F ,»e2o6si0e Leow sil Tin SS OFGY . GS OF 6 18° 
- ° _ 
7 np a Loe ’ n z ; ot a - 7 = 
bes Voters sideadit od? aeeviae fz. . Eteto of tasine sails -»s 
> - a ‘a , - ‘ 
P i a << iT, a ei xe ~ ie wy 7 
pebivety at. ypusm@e vivesta sAg ereinl 30° ee J 
’ 7 =n 9 - 


re 
wited £ Nea io yilsved myles 


— 
- sh Pe 7 
‘ 


+s abtaatog 


_ 


4 + i iad 
* 


in 


mT @. es 


e 
—< 


a 


sphere collision models(*2). For high energy incident ions 
the kinetic energy transfer is inelastic which involves the 
energy transfer not only in nuclear motion but also in 
excitation and ionization of the electrons of the target 
atom. 

The secondary electron emission yield is independent of ion 
energy when the ion energy is less than 1 KeV and potential 
ejection predominates in this region (*3). At energies 
greater than 1 KeV, kinetic ejection is not negligible and 
the yield increases, initially linearly with the kinetic 
energy of the ion and later linearly with velocity, 
eventually approaching a maximum, before decreasing (*4,*5). 
It was found by Hagstrum (*6) that for low energy ions, the 
neutralization process occurs with a high probability before 
the kinetic process can occur. The following neutralization 
processes of ions near the metal surface were studied by 


various researchers (*7,*8,*9,*10). 


a Resonance neutralization 
b Resonance ionization 

c Auger de-excitation 

d Auger neutralization 


The resonance neutralization was first suggested by Oliphant 
and Moon (*11) and the process is written as, 


“ - * 5 
A + N Em --------- A + ( N-1) Em -----(2-1) 


~ 
A 3; incident ion 


N ; the number of metal electrons 


7 

y > 4 
into z nine ena 
+ ee anche fabio be 
“ md. ets nue 


Y= ¢ 


aa } 4 


= =e 

sotsen 286i iowa at 7 we 

- ia _ . 7 na 

socal — ontasts. oily. a9 portes 

‘ ih, ih & 
| : - x 


7 ‘ie 


~ “4 
side tout 2eis yaar S262 of 


2. 
re) 


7 


“a” 
o 
po 
a 
7 
mal 
Pe 
= 
i 
2 
a 
: 
i=) 
~ 
ir} 
me 
oo 
a in| 
> 
h~ 
4 
So 
> 


faktaiod boe oan ) wade eagf-ci yp aa not oon weds eoteas. (| 
: ; » "4 . “4, 
tougye +A .(6*) notoed 2ia% es scimohery aor i ate 


. “5 f 4 ; » - = 7 
bos eidetpiltren ton.ek fotiosets 2578 «Vou ° aadt T3*8a%p 
4 i i +0 - 
Sue « 47 wa wre ies A i¢ 
risoatd edt dtiv viteert IDdsitin: (esasetoanr Skai¢ sdf 
ttcoofey dtiv yitsegil tatel ons aos ag 20 YeoTeqe 


(2¢ 04) polessTseb saoted (wueixes 6 uaiddsotags yalsa> wave 7 


- ; 
> % ~~ 2 “A - yt . i s ; -_ 
sit ,esol yoress wol wa tedt (2%) merrzepet ya fava s Says 7 
v . 7 - : ‘t 
stcted yeifidsdonad dpid « dtiw a2Fu900 2269019 dGiress lo ireen 
a - - hed * — <— —_ ~ = * os a ie: _ _ 
roitssiivztosa privoliot af? .a0os0 aie: ise f4 ad 


vd fotfpte sase sostzus Let 


> or or. | uskvAalgtsa iuven sa00k % on bb 
av q ‘ 
. a ear 
4 sig ilo vd hetcepanue saTti aew. datine tlatnyen posi +42 . 
of ts Pm se 

-” : - oS _ 

25 fTazsitr zt cosoon- of: pant 0 oem xl 

x sei 7 -Y 


(?=h)-- 22> et (7 a } * & ene at 4 ‘ + | 
. > 7 
: n : : a <{ - 4 7 wre ; 
; ¥ r 


10 
cf 
A 3; excited atom 
The process (2-1) can occur whenever the condition 2 < Em, < 
Wa is satisfied. 


where 


e9 ; average work function of the metal Em ; 
electrons in the metal 
Em 3; the energy of metal electron relative to the 
Fermi level Wi 
Wa 3; outer work function of metal 
The excited atom A may reverse the resonance neutralization 
process causing resonance ionization. This process is 
written as, 
a = 
R ---------- > AR + Em -------- (2-2) 
Also the excited atom in equation (2-1) may undergo Auger 
de-excitation . The process is written as, 
A + N Em ------- > A+ E + ( N-1.) Em ----(2-3) 
This process can occur whenever the condition eI > eg is 
fullfilled. 
where 
eI ; the ionization energy of the incident 
particle. 
Whenever the condition E’ > e> is fullfilled, secondary 
electrons, FE , are emitted. Eis excitation energy. The 
maximum kinetic enegy of the secondary electron ejected in 
an Auger de-excitation process is given by, 
( Ee )max = E = e--=------ (2-4) 


The resonance neutralization is followed by the Auger 
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de-excitation process to completely neutralize the ion. In 
contrast to resonance neutralization, which is a two step 
process, Auger neutralization is a single step process of 
ion neutralization on the conductive surface. The Auger 
neutralization process is written as, 

Ampeg NeEUMs==ss2-Esse2 >A +E + ( N-2.) Em ---(2-5) 
Two electrons from the conduction band of the target 
material are involved. One electron is used to neutralize 
the incident ion directly to the ground state and the other 
electron is ejected by the released energy from the 
transition of the ion. The maximum kinetic energy of the 
secondary electron from the Auger neutralization is given 
byy 

( Ee )max = eI - 2€9 ateetetetestetetsted (2-6) 
It is lower than that of the Auger de-excitation process. It 
was found by Hagstrum (*12) that Auger neutralization is 
more probable than resonance neutralization at low incident 
ion energies. Hagstrum's (*13) experimental results show 
that the maximum energy of secondary electrons ejected from 
Mo substrates irradated by ,He stays below ( Ee )max for 
Auger neutralization when the incident ion energy is below 
40 eV and for greater than 100 eV, the energy distribution 
of secondary electron was over (Ee)max for Auger 
neutralization showing resonance neutralization had occured. 
As suggested by Hagstrum (*14), the condition eI > 2€9 must 
be fullfilled to make Auger neutralization possible. From 


the work function and ionization data of carbon, lead and 
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magnesium, it is clear that the neutralization process of a 
singly charged low energy Pb ion depends mostly on the two 
step process, i.e., resonance neutralization and Auger 
de-excitation when the target surface is carbon or lead. 
When the target surface is Mg, the dominating neutralization 
process is Auger neutralization until the surface is covered 
by lead at which point the process changes into a two step 
one. For singly charged low energy Mg ions, the dominating 
neutralization process is a two step process when the target 
surface is carbon or lead but when the target surface is 
covered with Mg, the process turns to Auger neutralization. 
Table 2-1 (*15,*16) shows work function and ionization 


energy of carbon, magnesium and lead. 


2.23 Collison between the incident particle and a target 


atom. 


The ion beam energies employed for the experiment are 
less than 150 eV and the collison between a neutralized 
incident atom and a target atom is considered elastic. The 
hard sphere binary collison model can explain the kinetic 
energy transfer process. The total energy transferred to the 
target atom from the incident particle is given by (*17) , 

T=4M1M2/(M1 + M2)Eo sin@/2-------------- (2-7) 
where 
T ; The transferred energy 


Eo ; Incident particle kinetic energy 
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Table 2-1 


Work function and ionization energy 


ee ee ee Ss ee Se ee ee ee a ee ewe we ee 


ee 
eS 


work function (eV) 


w3 ew ew mM ww ew ew em = we we wm wow ww oe oe 


lonization energy (eV) 
Single 


double 


eee ee ee ee ee 


fexeumieyeysy Uf «vale; 
S=SSS====| —— = 
| 
4.39 een ae) 
R]— ee ewe ew @ | @2-ew eee @ w= = 
| 
ipi.ot 
| 
BPE T 
| 
{ 


om came um cee me ees ee ce ee wee 


13 


| 4 
nme Et erie pee ete amet Fey te me I 


ye, - 
ts me on 
' - 7 ! 7 
7 a Sess Pye az se 
Y¥Orojis 267 8h5 700% eae sods aie 7 
Oe ee ee ee, ee eed eerie oh eee ae 


- 


a a eae 2-2 oe ye 
A i 


PM o)-5 BOS Tes i Urs ; 

aii tae a= — SP2wSsees=e sce SeSse 

be Lane oe . . 
is ii 7 c j a Pr. E q ‘Yay $2On 243° dt ‘ j 

od ate Seems pale oe nl ere Se eee eee —— oa, ae & @ce 

ae i Way tes ane noksaxiaos 

pa,.T 7 b eofps £e 

SJ 

t 

“4 

4 


: 7 

. 4 ee v1 ' 7 ’ gj 7 
tee 2. ; é. > 

nay : } ; ( 

f' 


aoe 


' 

f 

: . tT 

Pau ib ms Po 
t 

i a 


———- <ee rw a 


Nes a4) = 


14 


M1, M2 3; Incident particle and target atom weight 

6 ; The deflection angle of incident particle 

after the collison in the center of mass system. 
When 6=180, iee., the motion of incident particle is 
reversed or stopped after collison, the maximum energy 
transfer is given by, 

Pe hth) Chia oT (2-8) 
In an elastic collison between the incident particle anda 
target atom, hard spheres of radius R are used to explain 
the process. R is determined by the interatomic potential 
and the energy of the incident particle. Several theories 
(*18,*19,*20) discuss the interatomic potential but none is 
definitely valid for all interaction distances. It is 
suggested that for small approach distances, a Bohr 
potential may be reasonably employed while larger separation 
distances require a Born-Mayer model. For example, an 
incident copper ion having a kinetic energy of 50 eV and 
colliding with a copper atom will have the following hard 
sphere diameters depending on the choice of theory. They are 
given by, 

0.722 A: Bohr model 

1.264 hs Born-Mayer model 


1.986 A: Thomas-Fermi model 


22.4 Nucleation and growth 
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LD 
The process of nucleation and growth of thin films was 
investigated theoretically by various researchers 
(*21,*22,*23) and the thermal-accomodation coefficient is 
defined as, 
ol, = (Ev - Er)/(Ev - E) = (fv - Tr)/(Tv - T) 
Soca 7 (259) 
where 
Ev ; incident kinetic energy of the incident atom 
Er ; energy of the desorbed atom before 
equilibration with the substrate 
E 3; energy of the desorbed atom after it has 
equilikrated with the substrate 
Tv,Tr,T 3; corresponding temperature 
It is suggested (*24) that complete thermal accomodation, 
ie. %=1, is obtained when the incident atom kinetic energy 
is less than the activaton energy for desorption(approx. 0.5 
ev). For thermally evaporated Mg and Pb atoms with kinetic 
energies of less than 0.17 eV per atom (equivalent to 2000 
re thermal accomodation coefficient is unity and all 
the impinging atoms are equilibrated with the substrate 
without any rebounced atoms. For ion beams with incident 
kinetic energy of 50 eV, it is very likely that the thermal 
accomodation coefficient is less than unity because the 
incident kinetic energy is about 100 times the energy 
necessary for desorption after equilibration with the 
substrate and this ratio is beyond the suggested maximum 


value of 25 (*25) . During the initial stages of film 
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deposition by an ion beam with a kinetic energy of 50 eV, a 


fraction of the incident atoms bounce off the substrate 


causing a lower initial adsorption rate. The adsorbed atoms 


can either stick permanently to the substrate or reevaporate 


in a finite time. The following stages of film growth will 


be discussed for the adsorbed atoms ; 


a 


The island stage 

The adsorbed atoms migrate over the surface of the 
substrate and they constantly aggregate and dissociate 
until the cluster grows to a critical size. Beyond the 
critical sized nuclei, the cluster does not dissociate 
but becomes stable and grows to a larger permanent 
island. 

The coalescence stage 

The islands continue to grow, and touch neighboring 
islands, coalescing to form a new island occupying an 
area smaller than the sum of the original two , thus 
exposing fresh substrate surface. New incident atoms 
are adsorbed on these freshly exposed areas and 
secondary nucleation occurs. 

The channel stage 

As the islands continue to grow and coalesce and before 
a continuous film is made, a stage is reached in which 
the film consists of islands a few microns in lateral 
extent separated by narrow exposed substrate material, 
iee.,channels, just a few hundred angstrom wide. 


The continuous film stage 
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The new nuclei formed in the channel area grow and join 
the large islands, bridging between the elongated 


islands and eventually filling the channels completely. 
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a The diffusion mechanism 
Diffusion usually involves one or more mechanisms 
depending on the structural defects of the material. 
The following are the mechanisms of diffusion. 
1 Vacancy mechanism 
This is the most common type of diffusion 
mechanism, In thermodynamic equilibriun a certain 
number of vacant lattice sites can be expected to 
be present in a crystal. Any atoms neighboring on 
a vacancy can then diffuse by jumping into the 
vacancy, the result being an interchange of 
position of an atom and the vacancy. 
2 Interstitial mechanism (Direct interstitial 
mechanism) 
An interstitial atom diffuses by moving directly 
from one interstitial site without causing net 
motion of any other atom. 
3 Interstitialcy mechanism (Indirect interstitial 


mechanism) 


‘eb ot wre yr “en 


are 


5 es a ; ” ea Mie boi 
| et ton ese ae Si: a ei 
— Retneands lit aware: — 
: nel, can fie’ a: : ; 
wievels ROD denis ode onkigex y ites 
ey soy OW 
eu | 7 Pave 4 ; 


, U Tan ; r 


: ; 7 ey 
sf . ew LEI o2Aasy goles ra 5 
[ f 2 5 Mae 
Y ’ : , y 7 e i] 


sat 


Th 6 eee . 
snetnedoom gs0m IO ano eevylovat 


agi aileie notentaib ed? a 


r 


Lsveu acta0ctzia. 


[niastsa 94% to atasted Estrzoeagarsty co onthaeqed 


roL2dIVFTLb ‘to sqvc adaie> %a00 eat aL, 22cF 

7 ~ 

aT ra? 

£523 ; meLadifl ives ofe seypomsed? «2. smabradoee 
r ‘on 

st betoon sd opD eerie es 7it3si 7EBOSY to T390a8 


aff offci prigeanst ya aeviPkb a: / AS Yor eQky 6) 
J eet be . : 
"So pesdaira rat "#8 eaee8 tig hel As. vt piers v 
— 


syouassy sd Sas wih a8 t 


febtitesetnd tois710) ue bostfoan ist tissizatad ig 

; : 7) , 

rr — a Se ainedoan | 
ld “ : ‘ ; i ' : : . : 

gicve:tb oulvom: yd eset last $itarota figs: r. 


t " rd ; ni, 


a 


_ ton palapao shoddy Fie rebtaab int = sont 
| - ¥ 0 ee 


' #8 7 i? 


, 029+ ante; Yas 9 to 
ae we Ae ale W he 
ARS easenaw' arn iia re 
io ree 7 PWS f at 
; | age ny | ae 


ia . 


ay 


t ; 


18 


The interstitial atom moves by pushing a normal 
lattice atom into an interstitial site and moving 
into the lattice site itself. The location of the 
interstitialcy moves twice as far as do either of 
the individual atoms. 

4 Exchange mechanism 
Two atoms in the lattice move simultaneously by 
exchanging places with one another. 

5 Divacancy mechanism 
Pairs of vacancies that are bound together as 
divacancies move causing diffusion. 

6 Crowdion mechanism 
When an interstitial atom is located at a site 
along a close packed direction the lattice defect, 
called a crowdion, is generated and an indirect 
interstitial diffusion mechanism can occur. 

a Relaxation diffusion mechanism 
Liquid-like disodered regions where atoms are not 
really bound to lattice sites may be generated in 
the neighborhood of a vacancy. Atoms located in 
this region can move through the crystal. 

8 Ring diffusion mechanism 
Three or more atoms jump simultaneously to the 


next site around a ring. 


2.6 Thermal diffusion 
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Because of thermal energy, the atoms in a metal crystal 
are in constant motion around the equilibrium lattice sites. 
Occasionally, an atom will jump to a neighboring site. The 
jump frequency increases as temperature increases. The 
probability of surmounting the energy barrier, which atoms 
must overcome to move to the next site, increases as 
temperature increases causing diffusion. A summary of 
diffusion theory based on the atomic jump will be presented 
in the following. The two dimensional diffusion equation is 
given by, 


te 
= G70 Kaa Vc we eew ome ome ee (2-10) 
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where 

J 3; the flux of the diffusing species 

D"; the tracer diffusion coefficient 

c 3; concentration of the diffusing species 

dc/dx ; concentration gradient 

y+ the atomic drift velocity from driving forces 

acting on the atoms 

The tracer diffusion coefficient, Di; and the atomic drift 
velocity V,is directly related to the jump frequencies of 
the atoms. The atomic concentration for the two lattice 
plane is assumed to be n1 and n2. The frequency of the jump 
from plane 1 is £12 and from plane 2 is £21. The total flux 


of diffusion is given by, 


eo ate en 2 ei mmm me (2-11) 
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where 
c ; concentration of the diffusing species 
a ; the distance between two lattice planes 
The local gradient in the concentration per unit area along 
the diffusion is given by, 
dod xeee (n 26-61) yanoc-acse-ss- = (2-793) 


From equation (2-10),(2-12) and (2-13), we can derive the 


following, 
D’ = a f 
‘Wo= a) (£12 - £21) <-=s=-=------ (2-14) 
where 
f = (£12 + £21 )/2 average jump frequency 


The tracer diffusion coefficient is proportional to the 
average jump frequency while the atomic drift velocity is 
proportional to the difference of the jump frequencies. 
Experimentally it is well known that the most important 
variables which affect Dare temperature, pressure and 
chemical composition. For practical purposes , the diffusion 
coefficient is written with the conditions of constant 
pressure and chemical composition as follows ; 

Dean exp( -E/kT) --+-------- (2-15) 

where 

Do ; experimentally determined constant for each 

diffusion couple 

E ; activation energy 

k 3; Boltzmann's constant 


T ; temperature 
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If the energy barrier height is not constant along the 
diffusion direction , the jump frequencies of the 
neighboring plane will not be the same and eventually the 
difference of the jump frequencies contributes to the atomic 
migration. The difference of the jump frequencies can be 
caused by the presence of one or more of the following ;: 
electric field, temperature gradient, non-ideal part of 
chemical potential gradient, centrifugal force from angular 


velocity or stress field. 
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CHAPTER 3 


LOW ENERGY ION BEAM DEPOSITION SYSTEM 


The low energy ion beam deposition system was designed 
and built by P. Bryce , J. Amano and R.P.W.e. Lawson in 1972 
for ion beam studies and since then slight modifications of 
the electronics system were carried out. The system 
specifications satisfy the basic requirements for high 
purity film deposition, which can be summurized as : 
a: low pressure in the substrate chamber; 
b : controllable deposition rate; 
c 3: readily adjustable final energy of arriving beam; 


d : minimum space charge expansion of the ion beam; 


3.1 Vacuum system 


The vacuum system consists of three interconnected 
chambers ; ion source chamber,ion beam chamber and substrate 
chamber. The normal operating pressure of each chamber is ; 
x 10 torr in the source chamber, 1 x 10 torr in the ion 
beam chamber and 1 x 10 torr in the substrate chamber. These 
pressure differences are realized through a differential 
pumping system. Figure 3-1 shows a schematic diagram of the 


differential pumping system. The calculated values for the 
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0.4 l/fsec, S2 = 100 l/sec, S3 = 900 l/Ysec, 
C1 = 0.1 1l/sec, C2 = 18 l/sec 
The pumps used in this system are as follows; 
For the ion beam chamber 
a 6 inch diffusion pump(S = 600 1/sec) 
b rotary pump 
For the substrate chamber 
a turbo molecular pump(S = 70 1/seéc) 
b Ti-ball sublimation pump(S > 1000 1/sec) 
c diffusion pump 
d rotary pump 
For the ion source chamber 
rotary pump ; used only for roughing out 
The partial pressure of residual gas in the substrate 
chamber is monitored by an A.E.I. Minimass residual gas 


analyser. The overall vacuum system is shown in Figure 3-2. 


3.2 Ion beam system 


The ion beam system is divided into three sections ; a) 
ion gun, b) ion beam transporting system, and c) decelerator 
and substrate holder. The ion gun had been modified from the 
basic Colutron system by J. Amano so that it can accommodate 


materials with a higher vaporization temperature. This is 
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31,$32,S3 ; Effective pumping speeds 


Figure 3-1 
Schematic diagram of differential pumping 
system 
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ION GUN SUBSTRATE | ,; 
CHAMBER CHAMBER . 
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Figure 3-2 
Schematic diagram of vacuum system 
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done by eliminating the problem of condensation between the 
furnace chamber and the discharge chamber. The schematic 
diagram of the ion source is shown in Figure 3-3. One of the 
major parts of the ion beam system is the ion beam transport 
system. The Colutron ion beam system Model G-2 was chosen 
for this application. This system consists of an 
acceleration and focusing system, vertical deflection plates 
and an E x B velocity filter. Figure 3-4 shows the complete 
beam system from ion source to target. 

The velocity filter consists of a magnet, deflection 
plates and guard rings. The E field is normal to both the 
ion beam trajectory and the B field to counterbalance the 
Lorentz force on a particle of given velocity. Under a 
balanced condition a selected velocity of ions will pass 
undeflected through the separator while particles with 
different velocities will be deflected. Figure 3-5 shows the 
geometry and mass separation in an E x B separator. In the 
uniform field region of the separator the motion of an ion 
will be described as follows when E and B are independent of 
SeOLaZ.« 

mx=eEx-e2 By 
mzz= eBy x 
BY eee (329) 
By solving the above equations for the balanced condition, 


1i.€. ,M x = 0, m y = 0, we get the special velocity Vo, 


The velocity of an ion entering the velocity filter is 
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Figure 3-5 
Geometry and Mass Separation in ExB Separator 
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determined by its mass, accelerating energy and ionized 


charge number. The velocity is given by, 


gq 
i 


(2n eV / MN) ----------- (3- 3) 
where 

M ; mass of ion 

V ; accelerating potential 

nh ; ionized charge number (n=1 for singly charged 

ion) 

e ; electronic charge 
The defined mass is selected when the special velocity Vo is 
adjusted to be equal to v in equation 3-3. During the 
experiment Ex is fixed to 5.0 x 10 Vm and only By is 
controlled. It resulted that the overall system mass 
resolution is M/gM = 4.5 and the detectable maximun 
molecular weight is 600. Also, the velocity filter is used 
to focus the ion beam into a line image instead of a spot in 
order to cover a wide area. Equation 3-1 describes only the 
ideal section of the velocity filter where the E and B field 
is ideal as shown in Figure 3-5. In fact, the E field is not 
ideally perpendicular to the Y direction but somewhat barrel 
shaped, because the E field is generated by a long parallel 
plate along the z-axis and the spacing between the plates is 
comparable to plate width. Any distortion caused by this 
design is corrected by the use of guard rings(Fig. 3-6). 
A desired ion beam entering the substrate chamber through 
the exit slit between ion beam chamber and subtrate chamber 


experiences first of all an electric field for horizontal 
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Figure 3-6 
Cross. Section of Velocity Filter 
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deflection. The horizontal deflection field and the 
decelerator( located slightly off the undeflected beam axis) 
cooperate to eliminate neutral particles that are created 
during the ion beam transportation process. Figure 3-7 shows 
this electrode structure. The decelerator plate is located 
53 mm from the exit slit and the exit slit is 110 mm from 
the velocity filter. The mass resolution for the velocity 
filter is given by (*2), 

M/AM =a lE /Y (4 VD) --------- (3-4) 

where 

a; the length of the velocity filter 

1 ; the distance from the center of the filter to 

target (drift distance) 

E 3; the electric field in V/n 

V 3; acceleration voltage or ion energy 

D ; the dispersion between mass M and (M- AM) 
Applying a=0.15 m, 1=250.5 mm, E=5.0 x 10 V/m, D=26 mm 
and we get the mass resolution, 

M/ofmM = 4.5 
With this level of mass resolution, the elimination of Mg 
isotopes from deposition is not likely to be possible but 
the detection of isotopes will become possible with a finely 
focused beam probe . A Van de Graaff accelerator which has 
He ion beam probe with 1mm diameter will enhance the mass 
resolution to 5 times and the detection of Mg isotopes 
becomes possible with this enhanced mass resolution. But the 


detection of Pb isotopes is not possible even with the 
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34 
enhanced mass resolution. As shown in figure 3-8, the block 
diagram of system electronics, the whole ion source chamber 
is floated above ground and kept at a high positive 
potential ( +4 kV). The accelerating electrodes are kept at 
ground potential. The decelerator plate potential is kept at 
the same voltage as that of the ion gun. An ion beam energy 
control power supply is inserted between decelerator plate 
and substrate with the substrate at the low potential by the 
amount of desired ion beam energy. Ion beam current on the 
substrate and decelerator plate are monitored separately and 
substrate current is integrated to control total charge of 


deposition. 
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Ton source anode power supply 
Acceleration power supply 

EFinzel lens power supply 

Ion source filament power supply 
Ion source furnace power supply 
Vertical deflection power supply 


w=Velocrty 11) ter BE field power supply 


Velocity filter magnet power supply 
Final deflection power supply 

Ton beam energy power supply 
Electron filament power supply 


Figure 3-8 
Block Diagram of System Electronics 
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CHAPTER 4 


RUTHERFORD BACKSCATTERING FOR THIN FILM ANALYSIS 


Since FE. Rutherford's nuclear scattering experiment was 
reported in 1911, an extensive range of particle scattering 
experiments has been studied by scientists. When an 
energetic particle beam irradiates a target, a fraction of 
the incident particles will be scattered back from the 
target by the Coulomb repulsion of the atomic nuclei 
(Rutherford scattering). The application of Rutherford 
backscattering to thin film analysis was begun in early 1970 
by Chu et. al. (*1,*2,*3) One outstanding feature of this 
technigue is the capability of depth resolution as well as 


mass, density and composition resolution. 
4,1 Mass resolution 


The mass resolution arises through the energy loss 
incurred by the small fraction of incident ions which are 
backscattered elastically by the Coulomb repulsion of atomic 
nuclei. The energy of the incident ion after backscattering 
is smaller than its initial energy Eo by an amount which 
depends both on the scattering angle and on the mass of the 


incident ion and a target atom. The relationship is given 
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by, 

E = Km Eo ------~--- (4-1) 

Km= (mcos6+ (M-msinOf)y/ (m+) - Sa — mean a (4-2) 

where 

E ; the energy of backscattered ion 

Eo ; the initial energy of incident ion 

Km ; kinematic recoil factor 

m 3; the mass of incident ion 

M ; the mass of target atom 

6 ; backscattering angle 
Knowing the parameters Eo, mand 6, and measuring the energy 
of the backscattered primary particle, the mass of an 


unknown target atom may be found. 
4,2 Depth resolution 


Most of the incident particles undergo small-angle 
deflections, penetrating deeper within the target. In the 
process of penetration, the incident particles lose energy, 
the backscattered energy of a particle re-emerging from 
underneath the surface after undergoing a large angle 
collision is smaller than that from the surface, because of 
its interactions with target atoms. The energy loss rate of 
ions in absorbing material is given by (*4 ), 


2 2 2 2 4 2 
-dE/dx = N Z (log2mv/I - log(1-B ) - B) ze /(4tEo 
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m,e ;electronic mass and charge 
ze,VvV ; particle charge and velocity 
Z 3; the atomic number of absorber 

N 3; the number of atoms jou 


E the average ionization potential of the ions 


ee 


BSs thesjratrvo w/c 
the negative sign in front of dE/dx signifies an 
energy loss. 


From a trigonometric analysis of Figure 4-1 ; 


E1 = Km Eo 

E2 = Eo - dE/dx t 

E3 = Kn E2 

E4 =9b3) Yd E/dx .t/cosd 


= Km(Eo-dE/dx t)- dE/dx t/cosé 

= Km Eo -dE/dx t(Km + 1/cos68) 
Therefore the energy loss 3 
AE = E1- E4 = (G@E/dx Km + dE/dx 1/cos®@ )t =(S) t 

where 

(S) = Km dE/dx)in + 1/cos6 dE/dx)out -------- (4-4) 
(S) is the backscattering energy loss parameter(*5). If we 
Ger heme (5) lOre ee ileands (S)5sbor 2 Ea) the ratio of (S), to 
(S), is given by, 
(S), /(S), a eee aces UE Cl ke) EmCEsck. Ti costck/dx) 


=(Km +1/cos6€,/ €,)/(Km +1/cos6 €, /E,) 


where €, ,€, and €, are stopping power at energies E,,£, 


and E, from Figure 4-1. 


’ 
en ha Poe to 


7 La 
7 a ee a ; : 
| oes +i ninok ‘Ss Giz: 2 
7 . atod sit 30 Sens at0q: toi ides intl > 966% veal : 3 
: 7 = et ; | oe ai oa eat a ; 
as eeibiipie BBY to tucat nf ania aby vay cae 7 
’ a ; My > cae vin ane a 
: _— $5 : 
x : t-! i ie ole sys 7 
ra i: 
‘Geonks cB\SH -— (2 “‘xB\S bot joe =. 
(Be0O\r + i ad #B\SD- a 
- ' ¥. 
: ener gerade od? osoteaud? 
> i 7 
+ 02) ae t nasa dl Pres Sb), = Ae - = : 
ve. owe. 
i ¥ ; i 
; eed aaecital PHO(KONTD 5 PeouNr fied a5: am — a 
. . mo a ad se 7 
a pr) 442 chaenanazag enol eae anneal feleiosh edt |f 
ee Pa a 
‘se Fae ‘er. ce 2 ie oh seme g 798 te hia "pe aot) = 
= a — - 
: 7 mi 7 Pa : 9a seas 
; ae . : - 
fi (ep doen nove LENE be az (OD ryshny 
oe ie an 
naga ahhh AN 2 NBR 
gh: ab 7 ine A fe hea ~ 
_" x. ws Laing revo ae : 
; oe e si o@ 9 as 


ene 
nae 


va = rte 7 : ; ‘ 7 ates ‘s 


) eo 
a ea 


OY 
For Mg film, the ratio is calculated by applying the 
following conditions: 

8=170 , t=5000A, Kmg=0.5171, E.=2.0 Mev,E, =1.0342 

MeV,E,=0.984 MeV, €,=45.11 x10 ,G,=55.52 x10 and & 

=55.36 x10 

(S);/(S), = 0.995 
For Pb film, the given conditions and derived ratios are : 

@=170 ,Kpb=0.9262,E,=2.0 MeV,,=125 x10,E, =1.8524 

x10MeV,€, =127.7 x10 , 

when t=1000A, 

E,;=1.8142 MeV, €, =128.65 x10 

(S), /(S), =1.022 

when +=2000A, 

E,;=1.776 MeV, €,=129. 14 x10. 

(S).6o)tesee 1055 
From the above calculations, it is concluded that (S) can be 
treated as independent of film thickness with less than 2.5% 
deviation when the film thickness is less than 50004 for Mg 
and 1000A for Pb. The above results are compatible with the 


reported deviation of 5% for a 5000A thickness in SiO (*6). 
4,3 Density resolution 


When a number of incident particles irradiate the 
target, only a small fraction of them are backscattered. The 
possibility of an incident particle being backscattered when 


passing through the target material depends on the number of 
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Figure 4-1 
Particle energy in the process of backscattering 
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target atoms per unit volume and the cross section of the 
target atom. From Figure 4-2, irradiating the target of area 
a and thickness dx, the total number of backscattered 
particles N is given by(*7), 
Ne e0/Aze) Om (Spe) dna ooo ee (6) 
where 
Q ; the total irradiation dosage of the incident 
particle in Coulomb 
z ; the number of ionization of the incident 
particle 
e ; electronic charge 


Q/{ze) 3; the total number of irradiating particles 


Zl 


; Avogadro's number 

p ; density of target material 

A 3; atomic weight of target material 

0 ; total cross section of target atom 
From equation 4-6, it is clear that if the density of the 
target material changes, then the backscattering yield will 
be changed. In practical cases the backscattering yield 
measurement is done only for a small solid angle, hence the 
yield expressed by a differential cross-section, rather than 
a total cross-section, is more useful. Equation 4-6 then 
becomes ; 

Ne -e@7(zcjed0/da (Noe 7A) dx do3-—-=2->—— (4-7) 
And the differential cross-section given by the Rutherford 
formula is (*7), 


a0/ao.= 1/4 (221 227(M1 -VIT_+)) 17sin (8/2) 


t4 j 


an ta 


ed? Yo coltaee eaoss ails ies omitov vinn zag aaptR sepa 


a sit 7 iia we : > 7 " : > eae : a 2 
989 3a tepist ‘ = = * ' - L4n4 — +i Tap? - 


nt 
Seresteae) 
LS 


‘Se p iene fs s Pe ota 2 7 Sea Miaka: TS 

; . 7 ig : 3 
yA } = ; | . (v yy y > a 2 2s faiix i 
7 p a 


——_ 


arene xin), af atone 
iu — + = _ 
ih | i we 
gashiogt edt to spscob nolksslsextt etre a navel 


ra in 
ee a 


duolvod ni sf trae 


tnebionlt ef? Fo noits ‘sgenet atti att ae 


: | | sfo £7289 


$p7eis sivorztsvels ; s 
> 7 - - - 


eelortiag paltaifenrii Jo isdapa Letor eit ; {eeyX\O 


wadlerin at ahs pews - 3 : 
: i * ye 
; LELTS IPE SETAE ¥e yi tensh 2 Oo ; 
feiastam +3022. 20 tipiow oigess 3 & 
wete teprst to moigoeae ebogn Lato? oO . _ 
: co i ’ 
ait to yitanceb off 2f teat — a£d. ee me yin to bkevpe aoa 7 
- : 
Litw Blety paites ssxeAS sd ent $ rods vasvanto Sefamias te sprst ' 
al 


biety paizestesadosd edt wai Reet a: me ea whee: weiD ad 


% a 
od 2 som: oi elpos Bifoe fflsae 8 r0T Yfae efoh at E stonsrvaon 


asd+ 196d362  tolkt>se-840 725 IsEa RLS ® yd hoes ota ) ek 


gad? a-O sottenpa .fvtedca ae ai aos tooaentor sii 
ie I Z i; me a 


a 


) | MN a _ 
ined) Hip Heat ess ‘6 AN Qi BENDE aa = 


7 a ; ey A ; ; [—< e ” 
: Ht034 editing = : ‘ e bie a0 eae aes orem s] & by . . S 
ye a #0, ; . ; ; ie ; 
- ia 7 : er mn 4 a a H ap 7 5) Te 


- BY (ks he ele Les 
orient! he mys ae ‘o) ON? = neue 


ery k 
ra i" piv ; ‘ip ape we. 


42 


Figure 4-2 
Backscattering yield from a slab of target material 
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where 
Z1,M1,V¥1 ; the atomic number, mass and velocity of 
incident particle 
Z2 ; the atomic number of target atom 
e ; electronic charge 
6 ; backscattering angle 
Equation 4-8 is expressed conveniently in terms of energy in 


MeV, as follows; 


- 2 
dQ/do,= 1.296 x 10 (21 D2 AR5° ) 1/sin (6/2) (b/sr) 


From equation 4-7 and 4-9, the density of the target 
material can be determined by measuring the film thickness 
dx and the total number of backscattered particles at given 


Q, Qa. and Be 
4.4 Composition resolution 


Composition ratio can be determined directly from RBS 
spectra. The total counts under each peak are directly 
related to equation 4-7. The ratio of the areas under two 
peaks is as follows ; 

N1 / N2 = (401 /dao Np1 /A1)/ (402 /da Noe2 /A2) 
The ratio of the number of atoms in a unit volume is, 

(Njo1 /A1) / (No2 /A2) 

=(N1/(d01 /da )) 7 (N2/(A02 sda )) 


=(N1 /N2) (22/21) -r--------------------- (4-11) 
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N1,N2 is the total count under each peak. From 4-11, 
measuring the area under each peak from the spectra, the 
ratio of the number of atoms is defined. In practical cases, 
the spectra is divided into certain energy steps. The 
backscattered particles that fall in an energy step are 
recorded as the number of counts in that channel. Then the 
number of counts in one channel is, 

Hm = Q/(ze) dO sda (Np /A)da 4x 

= Q/(ze) d0/da (Ne /A) da AE/(S) ------ (4-12) 
where (S) is backscattering energy loss parameter given in 
equation 4-4 and E is the energy width of the channel. In a 
Similar way, the ratio of the number of atoms in a compound 
is defined by the ratio of the number of counts in the peak 
channel as follows, 

(Np 1 /A1)/(Np2 /A2) = (Hm1/Hm2) (22/21) ( AE2/ 

AE1) --------- (4~13) 

The energy widths of each channel are all the same 

Se ARIS= 7 EZ 

therefore, Ni/N2 = (Hm1/Hm2) (22/21) -e2r----- (4-14) 
N1,N2 is the number of atoms of composition 1 and 2 in unit 
volume. Equations 4-11 and 4-14 are valid only if 
composition is homogeneous along the incident particle 
passage. If the composition is inhomogeneous and the 
composition ratio is changing from one dominant element to 
another without clear boundary regions( like the interface 
cf two films diffused into each other) special consideration 


for the calculation of the composition ratio is required. 
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4,5 Composition analysis of multilayer thin filn 


It is very possible that the interface regions of 
multilayer thin films consist of several elements making a 
transition region rather than a clear boundary between two 
films. The transition region may range from several atomic 
layers to several hundreds of atomic layers depending on the 
film production method, the treatment after and,or during 
film production and the film material combination. Within 
the interface region it is very likely that the composition 
ratio varies along the depth from one dominant element to 
another. In the process of the application of RBS thin film 
analysis techniques to the interface region of a multilayer 
film, the problem of stopping power and composition ratio 
calculation is considered. By measuring the film thickness 
of each of the constituents in a multilayer film, the 
overall average energy loss parameter (S) can be calculated 
from the RBS spectra. When a certain depth range within the 
film contains more than one element and the composition 
ratio has a gradient along the depth axis, the procedure of 
finding (S) by measuring thickness is not applicable because 
the stopping power dE/dx in the compound region is varying 
all along the depth axis. The stopping power for the 
multi-element region can be most precisely determined by 
measuring the energy loss and the film thickness but this is 
effective only when the composition ratio is constant 


everywhere in that region. Measuring the composition ratio 


6 pare ise sin ssois : 
7 to ibe 
owt denik attienl nia 8 raids: To 
visors. resange aos? pao ¥ ie nukpet i ok: 
ed? 20 prebaseeh aisyst oieors to ani 
onizvS 20,506 Jartse. taegicer? sar 


atdsiw .aottagtiso, i citotem afl 
4 ae | a A 


roltieognon sit) ted? yiotil yasv sé oa otpet eoeizerat. 
ie i 4 - _ 


” on 


at J ‘ - - pe - — ‘ 
of tnuemels tasataoh sao. Sos dtish. sa? caoke 2ei1sv © 
mali? sidd 248 Jo gotisotiqgs ed? Io-ses001g 967? AF »tenze 
. ; = _ : me) 
zavsit?tive a Yo agivpert enetietct edt ot esupindoet ekayls 
— , cy 
oitsi aoltimoqgoo bus yowog pRiqgote to e@eseldoxy oi> \oLt? 
_ F , 


: ; ne 2 me re 
zeegAntdat alid siz eatrone: = sbete shieaoo zt nottsiasls 


a a 

odt ,Mft2 szeyelssivas nae etumrrigeaas aA io doses 30 
"a u 4 

Lt wi 


bersinaiso ad aad (2) tetsessag @ 1919 sestovs ttsteve 


metoaga 808 o> ‘ena ; 
{ it vere 
er H> ont Sent becy’ 


to oruhenong 6 i+ 42 see adtqeb >t 
— 


' 
eanszed | odes foqe fon a pans ais secrets mae 


“a om 


ai 7 he -_ 
pabyssv. at woz nat nema ede ie - teu 
iy 7 


Sa/088 
" : ' r beat" ae sts 
ny ed? 2 2 #24) pagers: ads « : 


ne , tts 


f — 
t# sas tacii yaerkonag 3 reo oy 6 : woos 
= _ _ a 7 
a; x 
at: ake a aegoendsid 
th ; ; ioe? 


dante seit ‘ + ene, 
36% eons) wae > ce oA 


6 


in a certain depth increment, the stopping power of each 
increment slab can be redefined. Redefining dE/dx and 
composition ratio for each depth increment can be considered 
aS a probable method for determining the composition ratio 
profile along the depth axis. The whole process of 
redefining and recalculating dE/dx and composition ratio is 
done by the help of a computer. It is known that the only 
way to determine the energy loss parameter (S) when stopping 
power dE/dx is not available, is to measure the energy loss 
from the spectrum and the physical film thickness t+ (*8). 
The determination of (S) for an arbitrarily combined 
compound of 3 or 4 elements by the above method is not 
suitable for composition analysis because the measured (S) 
represents the overall average value of the measured 
thickness, which is usually measured from one side of the 
film to the other. This problem is somewhat alleviated by 
employing Bragg's additive rule of stopping power for a 
compound (*9). Bragg and Kleeman first postulated the linear 
additivity of atomic stopping power in 1905. Equation 4-14 
shows Bragg's additive rule written for a compound made up 
of three elements x,y,z having a composition ratio n1,n2,n3. 
The new stopping power for the compound is, 


C(x; ¥, cee ude xthen2cynhen 36z 


where 
€(x,y,Z) : the molecular stopping power of the 


compound. 


Th to ti 25 tayo a6: +Aane tod 
eye x6) ab ade :Sehom Sn 2686 97 " 
BessB ates ad ania aes od a bre tose 
eu: iievaehs tie sana 


ort ni oe 
ee Rojenspo as sfodw ei? selae Aaa: 


obs: to to qaioo at 


ak. ob re2 ‘nobtkeoquoe: bas +e\ab veut 


f = - 
ylao oft ted? nwomt eb EF , Tl ia i 5 2n ae nae anob- 
Lay 7 : 


j | oe 
niqqo?te seiw (2) wetemeteq sa0l rez9m@: in papers ad 


2aof ypisde sdz osvanen” sottentows fan st ONES ‘zaveg 


: ae! 


(ep s easndnbde aE fat. tig sd% fais susssage pits n032 


paaldaon ylz isitidas a6 10% (a) To no ttnaiars+ soar 
7 


ton e+ bodtem evods sd4¢ vd etcemsia a sof Fo ceil 


(2) Seyvessq edd eevsoed SEAEL SAS. MULTAROGIOS zo sidastee 
" : o 
bstunsen ea¢ to avlev spatere’ Listero off aoe veo2get 
7 - , 4 
sit 20 ebée ono aort besueasn videwew' et hoidw vesentolay 
b Al ery 
¥i horsivelis sadvenoe eit waltiowg ite sa3dto any 03, Pigs 


on 7 
Tt rewoOY poinaase 20 spas tibbs 2 ‘pps 74 prtyolaas 
7 _ : 


tseatl edt beteleseog 22452 cenest 


et «(2p hamodmes 
a a a 7 a 
Ure golseupe .0Ch ob sewog pad 
y ; bee : _ ¢ a x 
i ~~ - _ - 
qu obhse tavGgne0 5 3 gran wee  ovirisis ~ winds 
~< ¢ : : -_ 


‘9: - _ 
j mah ¢ — 
s€aySa,fan otter wotet ijmoD 5s paived x ‘aK >t rsenta > onsae 
_ : 7 i 


_o¥ : j 4) 2 ae 
ak Sanogans eae F yewog aad 
: i 


 S3Ed + vasa 
- 
, | : Be Riad 8 
: aa 1. eee 
opt 7e a 


ads Lo art 


> 46 
a “eo - 
‘ 7 ney ; / 
y 
yy) 


Big 


Ex,€y,E€z : atomic stopping power of x,y,z element. 
The validity of this equation has been checked by scientists 
from backscattering experiments. Feng et. al. (*10, *11) 
showed that the energy loss in Au-Ag mixtures was able to be 
predicted within 1% accuracy by Bragg's rule and Au-Cu 
alloys within 2% accuracy, and they proposed Bragg's rule is 
valid for all metal alloys. They also checked the stopping 
power of solid oxides and concluded that the stopping power 
of solid oxygen is 6 - 22 % lower than that of gaseous 
oxygen. In the process of composition analysis in the 
following sections of this thesis , the stopping power for 
solid oxygen is considered to be 85% of that of gaseous 


oxygen as proposed by Ziegler et. al. (*12,*13). 


4.6 Composition analysis applied to Pb-Mg-0O-C combination 


The atomic concentration profile of a compound of four 
elements,i.e., Pb,Mg,0,C, is obtained through the following 
procedure. This procedure, however, is only valid when the 
spectrum peaks from each element are separated. The RBS 
spectrum from Pb-Mg-O-C by 2.0 MeV He , gives four, 
well-separated peaks and it gives one reason why the film 
materials of Pb, Mg, and the carbon substrate are chosen 
initially. The definitions of all terms used in this section 
are given below : 

+ ; depth in cm(t=0 defines the surface) 


4t 3; depth increment 
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(S) ; backscattering energy loss parameter in 
eV/fcm (see equation 4-5) 
dE/dx ; stopping power in eV/cm 
€ ; stopping power in eV-cn/aton 
€in ; the stopping power in eV-cm/aton LOU 
particle energy of Fin 
Eout ; the stopping power in eV-cm/atom for 
particle energy of Eout 
Fin,Eout ; E for inbound and outbound particle 
Hm 3; the number of counts in a channel 
dd/fda ; differential cross section in b/sr 
subscript 1-4 ; designates element Pb,Mg,0O and C 
respectively 
Z 3; atomic number 
Zz 3; the number of ionization 
Q ; the total irradiation dosage of the incident 
particle 
AE ; the energy width of one channel 
e 3; electronic charge 
da ; the solid angle of detector effective surface 
measured from target 
nN ; Avogadro's number 
N ; no. of atoms in a unit volume(atomic density) 
Km 3; kinematic recoil factor 
The detailed procedure will be discussed in a later section, 
a summary follows ; 


ae Assume dE/dx for the surface of the film. If the 
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major elements of the surface are known, assume 
dE/dx of one of the major elements. 

Calculate the composition ratio of the surface 
using dE/dx obtained in step a. The composition 
ratio will be determined independantly from the 
error of the dE/dx for the surface assumed in step 
ae 

Correct dE/dx for the surface using the 
composition ratio obtained in step b. 

Recalculate the composition ratio of the surface 
using dE/dx obtained in step c.. 

Recorrect dE/dx for the surface using the 
composition ratio obtained in step d. 

Increase depth by fte Find out the backscattered 
energies for the each element located at At under 
surface. 

Assume dE/dx at the depth ft is same as that of 
the surface. Calculate composition ratio at At 
Correct dE/dx at ft uSing the result of step f. 
Increase depth by ft, now t=2 ft. Go to precedure 
g and h for 2 ft. 

Continue to increase depth by At each time in the 


same way until the desired depth is reached. 


4.6.1 The calculation of composition ratio 
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The composition analysis from RBS data is discussed by 
Chu eteal. (*14,%*15). It is indicated that the composition 
ratio obtained from an RBS spectrum is relatively 
insensitive to the errors involved in the values of stopping 
power dE/dx used in the calculation (*16). From equation 
4-12, the atomic density of each element is given by 
N1= Hm1 (S)1 K /(d01/da) , N2 = Hm2 (S) 2 
K/(d02/do) 
N3 = Hm3 (S)3 K/(d03/7do) , N4 = Hm4 (S)4 
K/(d04/da) ---------- (4-15) 
K = ze/(Q AE do ) 
The ratio of N1 through N4 is, 


NizN2:N3:N4 


Hm1 (S)1/(d01/do) : Hm2 (S)2(d02/do) 


: Hm3 (S)3/(d03/da) : Hm4 (S)4/(d04/daq) 


The ratio of differential cross section of each element at 
the same incident energy is given by equation 4-9 and 
equation 4-16 becomes, 


N12:N2:2:N3:2N4 


Hm14 (S)1/Z1 : Hm2 (S) 2/22 

| hcl eyo er 2 any (es) (eat eee ee ee ee (4-17) 
(S) is defined by equation 4-4 ,. It cannot be determined 
unless the ratio of compound element is found because the 
term dE/dx of equation 4-5 is indeterminate. Equation 4-5 
will be rewritten as, 


(S) = Km N €in + 1/cos® N €out 
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= Gin N(Km + 1/cos0 €out/Ein) ---2---e2---H-- (4-18) 
Figure 4<3 shows the backscattered energy from different 
scattering centers in a compound target. The energy loss 
parameter for each element in the compound is defined from 
this figure as follows : 

(Ss) = (215 - 1) /t 


(S)2 = (22. - B2)7t 


i 


(S) 3 (E3 =-E3) /t 
4 
(S)4 = (E4 - E4)/st 
Equation 4-18 is extended to a compound of four elements and 


the energy loss parameter for each element in the compound 


is given by ; 


(S)1 = €in N (Km1 + 1/cos@ € lout/Ein) 
(S)2 = €in N (Km2 + 2/cos8 E2out/Ein) 
(S)3 = €in N (Km3 + 1/cos@ €30ut/Ein) 
(S)4 = Gin N (Km4 + 1/cos® E4out/Ein) 
eorce---- (4-19) 

where 


Cin : the molecular stopping power of the compound 

for the particle energy of Eo 

Cout : the molecular stopping power of the 

compound for the particle energy of Kmn Eo 

N : the molecular density of the compound 

Kmn : the kinematic recoil factor for element n 
The ratio of (S)1 and (S)2 is written as, 

(S) 1/(S)2 


= (Km1 + 1/cos®@ Eloutyeéin) / (Km2 + 1/cos® €20ut/éEin) 
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Figure 4-3 
Backscattered particle energy from different 
scattering centers in compound target. 
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From equation 4-20, it is understandable that the ratio of 
(S) is insensitive to the systematic error of stopping power 
€& Figure 4-4 shows the stopping power € vs. energy (*17) 
for the four elements considered in the experiment. When 2.0 
Mev ,He beam is irradiating 5000 A thick target, the energy 
range of He” and the stopping power range of the target 
element is also shown and Fin, Eout represent the energy of 
4He in the path of ingoing and outgoing in the target. The 
energy range of oxygen is calculated by assuming the atomic 
density of oxygen is the same as that of the molecular 
density of MgO. The term €out/Ein in equation 4-19 is 
derivable from Figure 4-4 if the film is not composite. If 
the film material is composite, a new € vs. energy curve can 
be drawn for a certain composition ratio of the compound by 
the help of Bragg's additivity rule and €out/Ein from the 


composite will be found from the new curve. 


4.6.2 Depth profile of composition ratio 


The stopping power &€ for the surface of the compound is 
discussed in the previous section . The backscattering 
energy loss parameter (S) can be obtained if the atomic 
density, N, is given. But here we have another diffculty in 
handling a compound and it is to determine the atomic 
density N at each depth increment. The atomic density for 


pure bulk standard material of Pb, Mg, O and C, is given by 
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Table 4-1. Measuring the physical density of a compound, the 
atomic density can be determined but it is almost impossible 
to measure the physical density of the arbitrarily combined 
compounds in the depth scale. The depth profile analysis of 
an RBS spectrum is done by assuming that N for a pure 
element with standard density is still applicable when the 
element contains minor impurities in it. The error involved 
in the assumption of N will cause an error in the depth 
scale because the stopping power dE/dx (eV/cm) is given by, 

dE/dx = NG -----------------2-------=- (4-21) 
But the concentration ratio will not be affected by this 
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Table 4-1 


Atomic density for standard bulk material 
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EXPERIMENTAL RESULTS AND DISCUSSIONS 


Seo Production voftt thin’ film 


Both the low energy ion beam deposition and thermal 
evaporation techniques are employed to produce thin films of 
Mg and Pb on carbon substrates. During the process of low 
energy ion beam deposition, the target chamber pressure is 
Maintained below 1 x 10 torr to avoid possible impurity 
involvement during film growth. The ion beam is mass and 
charge filtered to exclude ionized impurities generated by 
the ion source, and the neutralized or multiply ionized 
particles from being deposited on the substrate. The 
selected ion beam energies for film deposition are 24, 48, 
72, 121eV for Pb and 24, 48, 100, 150eV for Mg. The 
combination of these energies makes 32 samples of double 
layer films of Mg and Pb on carbon substrates. From the RBS 
spectra of these samples, the most efficient film growth is 
observed at an energy combination of 48eV for Mg and 48eV 
for Pb ion beams in both film structures of carbon-Mg-Pb and 
carbon-Pb-Mg. (Under constant ion beam deposition dosage of 
0.22 Coulomb, maximum film thickness was obtained at the 


above energy combination) The optimum deposition energy for 
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double layer film of Mg and Pb is confirmed to be the same 
as that for single layer films of Pb and Mg. An incident ion 
beam energy of 48 eV was found to be the optimum for 
producing thin films of Mg and Pb(*2). The investigation 
focused mainly on films made at the optimum deposition 
energy of 48 eV. In addition to the above 32 samples of low 
energy ion beam deposited double layer films, 4 samples of 
double layer films of Mg and Pb on carbon substrates were 
produced by a combination of thermal evaporation and low 
energy ion beam deposition with ion energy of 48eV. Also, 2 
samples of double layer Mg and Pb were thermally evaporated 
cnto carbon substrates. The carbon substrates were all 
polished mechanically before being installed in the target 
chamber. The thermally evaporated films of Pb and Mg on 
carbon substrates were made in a vacuum bell jar kept below 
Sax There during deposition. The shapes of the ion beam 
deposited films were controlled to produce rectangular 
deposits of 1/4 inch wide and 3/4 inch high. Another 
deposition was carried out on the top of one deposited filn 
so that the two deposition patterns formed a cross. The 
center of the cross constituted the double layer film while 
each end of the bar comprised a single layer film of each 
material. The total dosage of deposited ions was set to 0.22 
Coulomb and the expected thickness was 3700 A for Pb and 
2800 R for Mg provided that the deposition was made very 
uniformly over the 1/4 x 3/4 inch rectangular shape, no 


sputtering occurred and the density of a deposited film was 
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the same as that of bulk material. In practical cases, the 
deposition shape changed as the ion beam energy was changed 
because space charge expansion increased as the ion beam 
energy decreased causing the deposition pattern to broaden. 
The film thickness therefore varied from the center area 
where the thickness was the highest to the boundary area 
where the thickness was the lowest. It was reported that 
deposition was negligible at Pb ion beam energies higher 
than 200eV(*1). For Mg ion beams a self sputtering 
coefficient of unity was observed at an incident energy of 
approximately 500 eV. The threshold energy for 
self-sputtering of Pb and Mg was reported as being 
approximately 24 eV (*2). If the ion beam energy is below 
the optimum deposition energy, sputtering yield becomes 
smaller while space charge expansion becomes larger. If the 
ion beam energy is above the optimum deposition energy, the 
former becomes larger while the latter becomes smaller. In 
the both cases, the film growth rate is reduced. The maximum 
growth rate is achieved only at the optimum deposition 
energy, i.€.,at the energy which has balanced sputtering and 


space charge expansion. 


5.2 Van de Graaff accelerator preparation for Rutherford 


backscattering analysis 


A He beam energy of 2.0 MeV was chosen because that 


energy is well below the possible resonance energy of low 
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atomic number elements (oxygen and carbon) involved in the 
experiment. The maximum stopping power and maximum cross 
sections are obtained around this energy, giving the best 
depth resolution. Figure 5-1 shows the RBS data handling 
system. The RBS spectra collected during the analysis was 
stored on magnetic tape and transferred to SDS 940 
minicomputer for further study. Final depth profile analysis 
was carried out by computer program written for the 
University's Amdahl 470 V/6. For backscattered particle 
detection, an ORTEC silicon surface barrier detector was 
employed and its specifications are as follows(*3) 3; 

active area ; 50 mm 

alpha resolution at 5.5 MeV 3; 22 KeV 

sensitive depth ; 1000 micron 

noise resolution ; 18 KeV 

operating bias ; -250 V 
Considering the dependence of energy resolution on the root 
mean square of energy, the energy resolution at 2.0 MeV is 
13.3 KeV. The secondary electron emission coefficients from 
Pb and Mg for 2 0 MeV He irradiation were measured as 3 


& Pb 


3.93 (electrons/ion) 

& Mg = 2.25 (electrons/ion) 
Figure 5-2 shows the effect of suppressor bias on target 
current during 2.0 MeV He irradiation for each material. The 
energy distribution of the secondary electrons ranges from 0 
to about 20 eV, with 85% from Pb (and 78% from Mg) having 


energies below 5 eV. These findings are consistent with 
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Figure 5-1 
RBS Data Handling System 
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63 
those of kinetic emission of secondary electrons reported by 
Kaminsky(*4). The minimum target current is observed when 
the suppressor bias is between -20 and -50 V. As the 
suppressor bias is increased to about -100 V, the target 
current slightly increased and then remained constant beyond 
a suppressor bias of -100 V. It is postulated that the 
minimum target current obtained at a bias between -20 and 
-50 V could infer a maximum energy of secondary electrons 
lying within this range. Figure 5-3 shows the detector 
geometry in the target room. ORTEC silicon surface barrier 
detector is placed at the one end of the cylindrical 
colliuator toj2uctease the) directionality, of detection. The 
resulting solid angle of the detector is calculated as 3.25 
x 10 sr. Theoretical backscattering yields of He from Pb, Ng 
and C are calculated by equation 4-7 and 4-8. The 
experimental backscattering yields measured from RBS spectra 
of Pb, Mg and C are compared with the theoretical values. 
The ratio of experimental to theoretical is 0.99 for Pb, 
1-63 for Mg and 1.04 for C. This result explains that the 
combined error of He total dosage, detector setting, the 
applied stopping power and the energy reading from the 


spectrum is less than 5 %. 
5.3 Stopping power and energy straggling data 


One of the most important parameters in the analysis of 


RBS spectra is the correct stopping power data for each 


“sopa89 oad an se ia! ot - boson 
- ey a = 


st moe ’ 
guys 


a ae 
a 


wr 
Bie» 


3 


- baoyed. beatae bentswex aedt - a ny = 


ot » sag nianees t20q ai $I mo o0r-4 
Ai 


i 
. 


1oitpeteh to wrils snot t ene att $269 


0b it ovode. £248 98 se 


oostava noohite I8TS9 smood tepret exit aici i 


hetsivolso ef totasted oft to eiprs Hiloe we 


[ 
¥ 


te alay ae hake 
ae m 1 be 
a ri ; oi 
Pal ook : 


co 
ais 


esid 6 te bomtstdo 
. ssh 
visbaovee to yptece meakzan aennt tie poe 


| ee de  ahdgew pe hei 


otzthaifys eft to bite ial ail as Sovelq af qos 


4) 


30 ehfeiy pe tant teoen oH 4 {ssfesne ant a ya OF =z 
io Ca > 1 a he Y : SI 
- é o oe ore | : 
i? 6-0 Boe v2 oe “4 beta ver aba era 26 5 


‘Fas . 
oy - 
Os a 
‘ ks 
41% vat ser 
we 
a 
Ve 
fi 
4 -= 


i 8 _ ibd tas 


“Aoa) a, re : 


64 


\£IZ0wW09y TO19Ee10Eq S_¥ 
€-S eansTy 


eIS sO pax So" ta eTSUY PTTOS 1049070q 
ee AS POEL INS CAt oo timo too) = 


YOLVNITIOO 


JW OS VauV 


qOvVdUNS YOLOALaC ea ewes 


YOLOALAC 4 NNN 


LHOYVL 


65 


target element. The stopping power data adopted here are 
contained in the Atomic Data and Nuclear Data Tables of 
Ziegler and Chu(*5). The adopted units of stopping power 
were taken as & (eV-cm/aton) unless otherwise specified 
(e-.ge AE/dx (MeV/mm or KeV/A) ). The unit of € is independent 
of density of the target material and provide a convenient 
calculation of stopping power of the compound. Table 5-1 
shows the stopping power data for alpha particles of between 
0.4 MeV - 2.0 MeV energy in Pb, Mg, O and C. 

The fluctuation of energy loss when energetic particles 
penetrate into target material is called energy straggling. 
The energy straggling theory of Bohr is modified by Lindhard 
et. al. (*6) 3; Chu has calculated energy straggling for He 
ions in various media using these modified theories. Figure 
5-4 shows the energy straggling for Pb and Mg in the 
backscattering experiment with 2.0 MeV ,He (*7) . The ingoing 
and outgoing paths are assumed identical and the ingoing He 
energy is assumed to be 2.0 MeV. The outgoing He energy is 
assumed to be 1.8 MeV for Pb and 1.0 MeV for Mg because the 
approximated kinematic recoil factor for Pb is 0.9 ard for 
Mg is 0.5. The energy straggling for low mass thin film 
material with less than 5000 A thickness is not serious 


compared with the detector resolution(13.3 keV). 


5.4 Analysis of Rutherford backscattering spectra 


Rutherford backscattering experiments were carried out under 
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Table 5-1 
Stopping RDA data (, He ) (*5) aS 
i gio dg ee tele hgin! wllnn eat Sarg Shae a tad ( x 10 eV-cm/aton) 
be 
| energy | Pb | Mg | 0 he Ree | 
| ( Mev ) | { | l { 
| SSS S333 22S ]2= | a | SSSSSS====| iS SSeS] SS | =S=SSSsS====}| 
| 1 { | l f 
| 2.0 } ese) {| 45.11 { 35.84 (pee 2 a 7 { 
| | | { (30.46) | { 
(aaaamaanaia eee |Senesesess [ces Geciaman a l 
i { | | { l 
| 1.8 i} aetsies 7 | 46.85 bes Sh (aez 7-08 | 
| | | | (32> 22)5} 1 
jSeenoechermec |pocoasaese | Saseesers Jnoeessoscs |ooeeseeee 
| | | { | ! 
| 1.6 { 132.4 { 48.82 { 40.27 | PRS | 
| | | | (34. 23) { l 
[soos es [nae peso cme [Pei (ero oso l 
| | | | { | 
{ 1.4 (e135. he Sie SS) { 42.81 yp SES { 
i { | | (36. 39) | { 
(Gee soo |PoStESose Cocoa aes jj opaoseeee l 
| j | { | | { 
| 1 2 ih Uae {536 2:1 (te to. 29 1e34.27 { 
{ | | I (38.50) | i 
joo Seeeeace (Seaman Ipeeseomses Joos [Poo | 
| | | { | { 
| 120 {| 141.6 i Bie is {! 47.34 Sis hs) ( 
I | | l (40.24) | { 
| souganoreoeen |nsesers Seal ppapease eos lamamaaes Soni) 2 ENS l 
{ | | | | i 
| Bless: } 141.9 { 56.78 | 48.39 yp OSs | 
| i | i (41. 13) | { 
| Seagate ee Peeoserre eat cea [poe es jecceatee 
I | { | | 1 
| 0.6 { 138.1 [Se AS | 47.72 [eee G58 | 
{ | | I (40.46) | { 
| =< ooo ee [oP see [pees sees [pS ee Jroosmescn- l 
| | l { | { 
} 0.4 {| 126.4 { 56.04 | 44.34 [oP 33.32 | 
{ { { { (37.69) | i 
| | { { 


* the value in parenthesis is for oxides 
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the following machine conditions. 


Ion beam energy : 2.0 MeV .He 

secondary electron suppressor bias : -200 V 
Ion beam current : adjusted in the range of 1 x 10" 
stax 10ecnoe depending on the target material to 
avoid piling up of data. normally for Pb 1-5 nA, 
Nog@507=871002n2 

Total ion dosage : 10 - 60 micro Coulombs. 
Backscattering angle : 170 degree 

Target annealing and cooling : target holder is 
heated by an electric heating element imbedded in 
the target holder. The target chamber pressure is 
maintained in the range of idetore during 
annealing. An iron - constantan thermocouple is 
used to monitor the target temperature. Argon gas 
is used to cool down the target to room 


temperature after each annealing 


£.4.1 Background level elimination 


Initially the backscattering spectrum is recorded into 4096 


channels by a Honeywell 516 minicomputer and later the 


channel numbers are reduced to 1024 by averaging each four 


channels. 


Before the spectrum is used as input data for 


depth and atomic concentration analysis, the background 
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level should be removed otherwise the background signal 
would give false input information for the analysis. The 
background level removal is done by the following peak 
finding method (*8) : 

The initial spectrum is assumed to be Gaussian and the 
spectrum is smoothed by using a Gaussian weighting function. 
The Gaussian weighting function has K times larger standard 
deviation than the initial spectrum. The value of K is 
empirically determined and in this experiment K=6 is 
selected. A new spectrum is made by convoluting the initial 
spectrum with the weighting function and the background free 
spectrum is obtained by deducting the new spectrum from the 
initial one. Figure 5-5 shows the procedure of the 
background level removal. 

From Figure 5-5, 

A: initial spectrum with background 


B 


background 
S : new spectra obtained by weighting function 
D : background free spectra 


The new spectrum is given by, 


si (@Wi Ak ) /eWj , where k= i+tj 
Djz-e 


S32 


Die Shia 


Ai: ith data of initial spectrum 


Wi : jth data of weighting function , Wj] = w(j) = 
As ae 


Ow : the standard deviation of weighting function 
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s : the standard deviation of the initial spectrum 


Ow=Ks, (where K=6 is chosen in this experiment 


By adjusting the value of K, the elimination level 
can be adjusted and set to the predicted channel 


axis zero intercepting pcint. 
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Figure 5-5 
Background level removal 
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5.4.2 CM(T) and CM(48) 


RBS spectrum of magnesium 

Figure 5-6 shows a portion of a complete 
spectrum for magnesium from CM(T) and CM(48). The 
leading edge of CM(T) lies ahead of CM(48) by 
about 50 KeV and the half maximum point of the 
front edge of CM(T) is 28 KeV ahead of CM(48), 
This energy shift is explained by the isotopes of 
Mg. Table 5-2 shows the Mg isotopes and their 
abundance. The calculated value of the energy 
shift between mass 23.985 and 24.986 is 28.4 Kev 
,and between 23.985 and 25.983 is 55 KeV. The 
thermally evaporated Mg film consists of mass 
232.985, 24.986 and 25.983 whilst the low energy 
ion beam deposited Mg film is made of single mass 
23.985. 

The oxidation of CM(T) and CM(48) 

The complete RBS spectra of CM(T) and CM(4U8) 
from which figure 5-6 are extracted, show an 
oxygen peak in both cases. Mg is chemically very 
reactive with oxygen and unfortunately the ion 
beam deposited films from the ultra high vacuum 
chamber have to be transported at atmospheric 
pressure to the storage bell jar, or to the Van de 


Graaff for RBS analysis. The oxidation of Mg has 
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been investigated by many scientists and the 
measured oxidation rate of Mg is given in Table 
5-3 (*9,*10). If we assume a magnesium oxide layer 
is formed on the surface of the Mg film and the 
slab of MgO grows thicker as oxygen is taken up, 
the oxide layer growth rate can be calculated from 
the above data. The oxidation rate of Addiss(*9) 
is equivalent to 2 A 7/BINS.- 4.7 A /mine and 
Cohen's(*10) data is equivalent to 0.4 A RID. = 
Ze X /mMin. provided that the magnesium oxide has 
the standard value of density. According to the 
above data, the storage of pure Mg films in the 
vacuum bell jar which has about 10 - 10 LoL 
oxygen partial pressure, cannot keep the sample 
from oxidation beyond several days. The 5000 A 
thick Mg film takes only about 30 - 200 hours to 
be oxidized in the vacuum bell jar at room 
temperature. Experimental results of the oxidation 
rate of Mg (*11 ) show that the sticking 
probability of oxygen does not drop when thicker 
oxide is present. This finding is not expected on 
the basis of the usual concepts of oxide film 
formation. This phenomenon is explained by 
Landsberg (*12) who assumes that oxygen is 
chemisorbed on active sites, which are 
deactivated, but create at the same time new ones. 


Figure 5-7 and 5-8 show the derived depth profile 
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and atomic concentration of CM(T) and CM(48). 
Landsberg's assumption is very compatible with the 
oxygen and magnesium concentration profile which 
shows fresh Mg on the surface and no oxide surface 
barrier which would block the intake of oxygen. 
Oxygen is evenly distributed in the Mg region 
rather than forming a high oxide concentration 


region in the surface area. 
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Depth profile of CM(48) 
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5.4.3 CP(48)M(48) and CM (48) P (48) 

Figures 5-9 and 5-10 show typical RBS spectra 
and Figures 5-11, 5-12 show the derived depth 
profiles of the CM(48)P(48) and CP(48)M(48) film 
system. From both depth profiles, no sharp 
interfaces exist between adjacent material. The 
composite region between the Mg and the Pb layer 
ranges over 1000 Kh wide » The composite region 
between carbon and Mg or carbon and Pb layer is 
much narrower than that of the Mg-Pb interface. 
The penetration of Pb atoms into the Mg layer is 
more noticeable than that of Mg atoms into the Pb 
layer. The probability of preferential diffusion 
of Pb into Mg is seen from the surface topography 
of both films. Figure 5-13 shows SEM topography of 
CM(48) and CP(48). The surface of CP(48) shows 
grains which are much larger than those seen from 
the surface of CM(48). It is suggested by Shewmon 
(*16) that the major atomic migration route is via 
grain boundaries when the temperature is below 
0.75 T(melt) K. Campisano et.al. (*14) pointed out 
that the preferential diffusion direction was 
controlled as the grain size of the host film and 
receiving film was altered. The preferential 
diffusion of the film with larger grain size into 
the film with smaller grain size is observed and 


it is explained by a longer grainboundary path 
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existing in the film having a smaller grain size. 
The surface topographies of the films illustrate a 
considerable amount of holes and cracks on the 
film surface. These holes and cracks are checked 
by the SEM in the electron probe mode and the 
characteristic X-ray coming from these holes and 
cracks do not show any difference from those 
coming from smoother film areas. The holes and 
cracks are only on the surface and they do not 
reach to the bottom layer or substrate. The 
surface roughness of the film could be due to the 
surface roughness of the carbon substrate. Figure 
5-14 shows the surface topography of carbon 


substrate. 
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Depth profile of CM(48)P(48) 


30 40 
DEPTH(X100 ANGSTRGM) 


84 


50 


— 
=~ 


‘¢ 4 


a ae 
; ; 7 
| Mule aOR + i 
1% AIX: 2 | 
. 7 


a 
. 


‘a ay ao 


a) 7 . nme 
i ta WUSAHO ft 
O : 7 


s 

: ae —_ 
a 
x 


ap P 
i" = 

mi wed 

ar 
A = 
i _ 2 
} o a : ar 
i? if f “> 


* 
-~ 


4 
ey 
—_ 


{ 
th 


= 
ha! : 
° a a ry 
. = i Ff 1 4 bet , 
iN J i = 
/ —_ 
D 
' aay 
. ‘ \ ” ‘ 
i ih Pa ‘= 
a - ” 
. ” ed 
] p ; 
U 
win. 
J oe a : 


= > 

a’ 
' ~ eel 

- ae 


( fo, # " Fa ea 
| pv, 1 Ky 
: z : ‘ 
j ‘ 7 ‘y \ 
eg -_ ‘ou 
: aS. 9 ’ ) Vi 
eee {= SE ees ee ee 
On 4 
‘Oe LP TS os ; GF. - , 
7 - & ¢ re mo “if py e ‘ ¢ aa ; ” 
GAG ¥ b Ou viet Lt} [xX a) ; 130 : ‘ : ¥ “308 "Wya iy 


; <aee 


i? 


of (34 )2(64)Ma 30; af ign 
ee a } 


ial 
7 i 
ae 

i) 


ea 
Ty 


. i ai . : ; 


ATOMIC CONCENTRATION IN Z% 


100 


80 


60 


40 


20 


85 


CP(48)M( 48) 


LEAD 
MAGNESIUM 
OXYGEN 
CARBON 


0 10 20 30 40 50 
SURFACE DEPTH(X100 ANGSTROM) 


Figure 5-12 
Depth profile of CP(48)mM(48) 


a > 


— - _ 
, 4 < Eee — 
TRAnTaaS 2 SInOTR 7: . — 
7 _ = . ‘. 7 7 = -_ 
a ; ee  -f- 
7 : ~ 2 meats a 


. 

| 

0 
beige 


4 
, 
LF | 

} 
| 
| 
J 
ia 
| 
yo 
om 
ae 
ees) 


4 
j 
ATOM OOTY 4) 


Oy 

d 

i 

+ 
aT 

i] é 


i 


By) oe. cuit (96 Woe bo ‘yridetgoret 


ee 
wa 


| oF 
1 
i 
. 
rc 
SB 


Figure 5-13 
SEM surface topography of CM(48) and CP(48) 
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Figure 5-14 
EM surface topography of carbon substrate 
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5e4.4 The surface topography of thermally 
evaporated and low energy ion beam deposited 


films. 


Figure 5-15 and 5-16 show the surface topography 
of CM(T) and CM(48), and CM(T)P(T) and 

CM(48) P(48). A typical distinguishing feature 
between thermally evaporated films and low energy 
ion beam deposited films is the grain size. The 
surface of CM(48)P(48) shows grains of about 2000 
A in size but that of CM (T)P(T) shows no 
distinguishable grains. Between CM(T) and CM(48) 
the distinction is seen in the same way, ieee, 
CM(48) shows very small size of grains but CM(T) 
shows smoother surface without any observable 
grains on it. No proof is advanced to explain the 
differences in the surface structures as a fuction 
of deposition method but it is postulated that the 
differences in incident particle energy, growth 
rate and substrate temperature during deposition 


by the two methods are the important parameters. 
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Figure 5-15 
SEM surface topography of CM(T) and CM(48) 
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Figure 5-16 
SEM surface topography of CM(T)P(T) and CM(48)P(48) 
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52.4.5 Annealing of CP (T)M(48) 


CP(T)M(48) was annealed with four annealing 
steps(for 30 min. at 100 C, 20 ‘minapats2001GA030 
mins at 300 C and 250 min. at 320°C). After each 
annealing step, the target is cooled down to room 
temperature by injecting Ar gas into the target 
chamber. A RBS spectrum was taken after the target 
temperature had been stabilized to room 
temperature. Figure 5-17 shows the RBS spectra at 
room temperature, after 30 min. annealing at 300 c 
and after 250 min. annealing at 320 ¢. Figures 
from 5-20 to 5-24 show the derived depth profile 
for the as deposited and the four annealed ones 
listed above. A small quantity of lead is observed 
to move into Mg region after 100 C - 30 min. and 
200 - 30 min. annealing step. A greater diffusion 
of lead is observed after the 300 Cc - 30 min. 
annealing step. After the 320 Cc - 250 min. 
annealing step, lead appeared on the surface of 
the film and this two layer region made a compound 
region of Mg and Pb whose composition ratio ranges 
front i? 1 to Sno ae diffusion of Mg into Pb 
region is noticed after 300 C and 320 C annealing 
steps but the diffused quantity of Mg into Pb is 
much less than that of Pb into Mg. From depth 


profile diffusion coefficient of Pb into Mg or Mg 
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into Pb can be calculated. For example, the number 
of Pb atoms moved into Mg during 300hce= "30 min. 
annealing period is derived by comparing the two 
depth profile taken before and after annealing. It 
is assumed that the diffusion is carried on only 
during the annealing period and the small amount 
of diffusion known to be ocurring even at room 
tempe rature(*15 ) is ignored in this calculation. 
Figure 5-25 shows the approximated concentration 
profile of Mg - Pb interface region before and 
after 300 C - 30 min. annealing. The derived 
diffusion coefficient of Pb into Mg at 300 c acy 
Whos Py ee) 5-4 TOM ena sec ) 

The detailed calculations are given as, 
The concentration gradients before and after 


annealing is, 


cad 
dc/dx = 1.806 x 10 atoms/cm’ a before 
annealing 

2) 
dc/dx = 1.09 x 10 atoms/cm eaten after annealing 


The total Pb atoms moved into Mg during annealing 
period is, 

dex ta oe ex te eri tae 

a co ay PAP ee 10 atoms/cm=sec 

Diffusion coefficient, D, is given by, 

De=) J 7) (dc/dx) 

taking the average value of dc/dx, 


27 
dc/dx = 1.448 x 10 atoms/cm* 
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Di= 10236 x folcuvees 
The preferrential diffusion of Pb into Mg is observed after 
200 C - 30 min. annealing step. No significant diffusion of 
oxygen is observed after each step of annealing and the low 
diffusion of Mg into Pb could be the cause of low diffusion 
of oxygen because the oxygen in the film is in the form of 


MgO. 
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Figure 5-20 
Depth profile of CP(T)M(48) : as deposited 
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Depth profile of CP(T)M(48) 
annealing at 100C. 
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Depth profile of CP(T)M(48) : after 30 min. 
annealing at 200C. 
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CP(T)M(48) AFTER 300C 30MIN. ANNEALING 
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Figure 5-23 
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annealing at 300C. 
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Depth profile of CP(T)M(48) : after 250 min. 
annealing at 320C. 
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Figure 5-25 
Pb atomic concentration slope in Mg region 
before and after 300C-30min. annealing. 
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CHAPTER 6 


CONCLUSIONS 


This chapter summarizes the present study of low energy 


ion beam deposited multilayer thin films. First the 


advantages of low energy ion beam deposition method compared 


to the conventional methods are summarized and then 


Rutherford backscattering analysis for compound material is 


discussed. 


Finally the investigated characteristics of the 


multilayer thin films are summarized. 


Ae | 


The advantages of low energy ion beam deposition 
method. 

The definite advantage of this method is that it 
has more controllable deposition parameters than 
any other conventional methods. The controllable 
deposition parameters are incident particle 
energy, film growth rate, particle directionality, 
pressure during film growth, mass selectivity and 
particle status. The controllability of deposition 
parameters gives the controllability of deposited 
film characteristics. 

Rutherford backscattering analysis for compound 
material 

Depth vs. composition ratio profiles are derived 
from RBS spectra of compound material of Mg,Pb,O 


and C. A computer program is written for the 
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derivation of physical profiles from RBS spectra 
having separated multiple peaks. 

The investigated film characteristics 

Both RBS spectra of thermally evaporated film and 
low energy ion beam deposited film show a high 
degree of oxygen contamination. This is considered 
to be a post deposition effect. The involved 
oxygen does not stay on the surface of the Mg 
films forming magnesium oxide layer but is 
dispersed into Mg region. 

RBS spectra from low energy ion beam deposited Mg 
film show that isotopes are removed. 

A preferrential diffusion of Pb into Mg is 
observed and the difference of grain size between 
Pb films and Mg films is considered to be the 
major cause of this result. 

The derived width of the interface region between 
Mg and Pb film layer ranges approximately from 
1000 A to 3000 rf and they are not dependent on the 


film deposition method. 
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